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Abstract. We have calculated the (sub-)mm spectral indices of 26 Herbig Ae/Be stars, for which we can determine 
the infrared spectral energy distribution (SED). We find a clear correlation between the strength of the ratio of 
the near- to mid-infrared excess of these sources, and the slope of the (sub-)mm energy distribution. Based on 
earlier multi-dimensional modeling of disks around Herbig Ae stars, we interpret this as a correlation between the 
geometry of the disk (flared or self-shadowed) and the size of the grains: self-shadowed disks have, on average, 
larger grains than their flared counterparts. These data suggest that the geometry of a young stellar disk evolves 
from flared to self-shadowed. 
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1. Introduction 

Herbig Ae/Be (HAEBE) stars are the somewhat more 
massive analogues of the T Tauri stars, which are low- 
mass young-stellar objects. The spectral energy distribu- 
tion (SED) of HAEBE stars is characterized by the pres- 
ence of an infrared (IR) flux excess, due to circumstellar 
dust and gas. The geometry of this circumstellar matter 
is believed to be disk-like (e.g. Mannings & Sargent 1997, 
2000; Fuente et al. 2003). 

Meeus et al. 2001 (henceforth MOl) classified 14 iso- 
lated HAEBE sample stars into two groups, based on the 
shape of the SED. Group I contains the sources in which 
a rising mid-IR (20-100 /im) flux excess is observed; these 
sources have an SED that can be fitted with a power- 
law and a black-body continuum. Group II sources have 
more modest mid-IR excesses; their IR SEDs can be recon- 
structed by a power-law only. MOl suggested phenomolog- 
ically that the difference between the two groups is related 
to the disk geometry. 

DuUemond 2002 (henceforth D02) and DuUemond & 
Dominik 2004 (henceforth DD04) have modelled young 
stellar disks with a self-consistent model based on 2-D 
radiative transfer coupled to the equation of vertical hy- 
drostatics. The model consists of a disk with an inner 
hole (~0.5 AU), a puffed-up inner rim and an outer part. 



The outer part of the disk can be flared (as in Chiang & 
Goldreich 1997), but can also lie entirely in the shade of 
the inner rim. The SEDs of fiared disks display a strong 
mid-IR flux excess, while self-shadowed disks have a much 
steeper slope towards long wavelengths. D02 can explain 
the difference in SED shape in HAEBEs (as expressed by 
the classification of MOl) as the result of a different disk 
geometry; flared-disks systems have group I SEDs, sources 
with self-shadowed disks can be linked to group II objects. 

In this paper we calculate the spectral index of the 
(sub-)mm (350-2700 /im) SED for a sample of 26 HAEBE 
stars. We look for a connection between this index and 
the geometry of the circumstellar disk. In Section |21 we 
classify the sample sources and present the computed 
(sub-)mm spectral indices. Interpretations of the results 
obtained in the latter section are given in Sectional We 
summarize our conclusions in the final Section^ 



2. Observations and Analysis 

The sources were selected based on Table 1 and 2 in the 
latest catalog of HAEBE stars (The et al. 1994). We en- 
larged our list of HAEBEs with objects from the sample 
studied by Malfait et al. 1998 that satisfy the criteria pos- 
tulated by The et al. We excluded sources that have the 
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Fig. 1. The SEDs of the sample sources with sufficient data to determine n. The circles indicate photometric mea- 
surements with error bars. The full line represents a reddened Kurucz model corresponding to the photosphere of the 
central star. The model was fitted to the optical photometric data. The (sub-)mm slope of the SED is given by the 
dashed line. 



amorphous 10 micron silicate feature in absorption (Acke 
& van den Ancker 2004; henceforth AV04), since these 
sources are probably in an earlier evolutionary stage and 
thus highly embedded. Early-B stars are excluded as well, 
because the disk-like nature of the circumstellar matter 
in these sources is still a subject of debate (e.g. sec re- 
view by Natta et al. 2000 on HAEBE disks in the book 
Protostars and Planets IV). For the remaining set of stars, 
we gathered broad-band J, H, K, L and M photometry 
and IRAS 12, 25 and 60 micron measurements in the lit- 
erature. For 26 of the stars for which this IR photometry 
was available, also (sub-)mm observations wc;rc; present in 
the literature. Our final sample consists of these sources. 

At (sub-)mm wavelengths, one observes the Rayleigh- 
Jeans part of the SED of the cold grains. Grains that are 
large compared to the wavelength at which they radiate, 
have a spectrum A^a cjc A~" with spectral index n = 3 
(black body). When the grains are smaller, n is larger, if 
the dust is optically thin at these wavelengths. Hence the 
(sub-)mm spectral index can be interpreted in terms of 



the grains size distribution of the cold dust in the circum- 
stellar disk (Sylvester et al. 1996, van den Ancker et al. 
2000, MOl). Nevertheless, this should be done with some 
caution, as we will discuss in the next section. We com- 
pute the (sub-)mm spectral index fitting a power law to 
the photometry between 350 and 2700 /im. In case only 
one (sub-)mm photometric point was available, we deter- 
mined a lower limit to the (sub-) mm spectral index us- 
ing the (sub-)mm data point and the IRAS 100 micron 
measurement. For BF Ori, we used the IRAS 60 micron 
measurement, since the IRAS 100 micron measurement 
resulted in an upper limit in this case. 

We note that the fluxes used to compute n come from 
a variety of literature sources, both single-dish and inter- 
ferometric. However, the vast majority of measurements 
come from single-dish instruments. There are only four 
sample sources (V376 Gas, EUas 3-1, HD 34282 and UX 
Ori) in which interferometric data are used. The latter 
two sources (HD 34282 and UX Ori) appear rather iso- 
lated (Malfait et al 1998). Therefore the interferometric 
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Fig. 2. Diagram based on van Boekel et al. 2003. 
The ratio Lnir/^ir is plotted versus the non color- 
corrected IRAS [12]-[60] color. The dashed line empiri- 
cally separates the group II (left) and the group I (right) 
sources. The size of the plot symbol is proportional to the 
(sub-)mm spectral index n. Open plot symbols indicate 
lower limits. 



fluxes are expected to be identical to those obtained with 
single-dish measurements. For the two sources with possi- 
ble envelope emission (V376 Gas and Elias 3-1; Henning 
et al. 1998, di Francesco et al. 1997) for which we have 
mixed interferometric and single-dish measurements, we 
do not note any systematic differences between the two. 
In Fig. n the SEDs of the sample sources with sufficient 
(sub-)mm photometry to determine n are presented. 

The IR SED of the sample sources can be character- 
ized with two quantities: the ratio of the near-IR lumi- 
nosity -Lnir (derived from the broad-band J, H, K, L 
and M photometry) and the mid-IR luminosity Lm (the 
corresponding quantity based on IRAS 12, 25 and 60 mi- 
cron measurements), and the non-color corrected IRAS 
[12]-[60] color. Both parameters compare the near-IR pho- 
tometric data to the mid-IR IRAS measurements (van 
Boekel et al. 2003; DuUemond et al. 2003). The near-IR 
SED is expected to be similar for all HAEBEs (Natta et 
al. 2001), while the major differences occur in the mid- 
IR SED. The luminosity ratio Lnir/-^^ir represents the 
strength of the ncar-IR compared to the mid-IR excess, 
which is lower for group I than group II sources. The shape 
of the mid-IR SED of a group I source is "double-peaked" 
compared to the SED of a group II member. The IRAS 
[12]- [60] color index provides a quantitative measure for 
this difference in SED shape. Group I sources are redder 
than their group II counterparts. 

Following van Boekel et al. 2003, we use these quan- 
tities to distinguish between group I and II in the clas- 
sification of MOl. We plotted all sample stars in Fig. [21 
The dashed line in the diagram represents iNm/^iR = 
([12] — [60]) + 0.9, which empirically provides the best sep- 
aration between the two groups. The sources defined by 
MOl as group I and group II do indeed occupy the correct 



location in this diagram using this separating line. The 
diagram provides a quantitative alternative to the subjec- 
tive classification method of MOl. In Fig.|21 the sizes of the 
plotting symbols are proportional to the (sub-)mm spec- 
tral index. 

The computed (sub-)mm spectral indices n of all sam- 
ple stars are listed in Table The error on the spectral 
index a„ reflects the dispersion on the fitted line. Note 
that the values of n show significant variations from source 
to source. The (sub-)mm spectral indices range from the 
default ISM value of optically thin dust n « 5 (Miyake 
& Nakagawa 1993) to the optically thin black-body value 
n = 3. The average values for group I and group II are also 
indicated in Table □ The error on the average is given by 
'^/'^n,j)~^ in which the an,j represent the errors on 
the individual n- values in the group. 

3. Discussion and Interpretations 

3.1. Comparison with other articles 

Hillenbrand et al. 1992 (henceforth II92) also classified 
HAEBE stars in groups. We stress that this was per- 
formed based on the near-IR (~3.5-10.6 /im) spectral in- 
dex. Their classification is not linked to the MOl classifi- 
cation used in this paper; both H92 group I and II contain 
MOl group I as well as group II sources. H92 group III is 
not relevant for sources with significant dust disks. The 
classification of MOl focusses on the ^--^20-100 /im SED, 
which makes it well-suited for this observational analysis. 

There exists a large variety of models for the SEDs 
of HAEBEs in the literature (e.g. H92; Natta et al. 1992, 
1993; Mannings 1994). While most of these (disk and/or 
envelope) models give satisfactory fits to the mid-IR SED, 
they have difficulty to accurately describe the near-IR 
bump that is typical for HAEBEs. One of the main im- 
provements in the model of DuUemond et al. (2001) and 
DD04 is the inclusion of the puffed-up inner rim, which 
generates this near-IR bump, as was first suggested by 
Natta et al. (2001) and Tuthill et al. (2001). Furthermore, 
the DuUemond et al. model explains in a natural way the 
roughly two types of IR SEDs (MOl group I and II) ob- 
served in HAEBE stars. The quantitative fine-tuned mod- 
elling of the SEDs of our sample stars, however, is beyond 
the scope of the present observational article. 

3.2. Grain growth in circumstellar disl<s around 
HAEBEs 

If the disks surrounding HAEBE stars are optically thin 
at (sub-)mm wavelengths, the spectral index of the SED 
in this wavelength range will be independent of disk ge- 
ometry. In that case, intrinsic properties of the cold grains 
in the outer parts of the disk must be responsible for the 
observed variations in the (sub-)mm spectral index. 

In order to test the hypothesis that the circumstel- 
lar disks in our sample are indeed optically thin at 
(sub-)mm wavelengths, we have computed the radii of 



Table 1. The (sub-)mm photometry and spectral index of our sample sources. The sources are sorted by group. The photometry gathered from the literature 

is listed in columns (3) to (10). The unit is Jy. Between brackets, the observing wavelength in micron is given. In column (11), the disk mass estimated based 
on the 1300 micron flux is presented. The (sub-)mm spectral index n and error c7„ are given in column (12) and (13) respectively. The references are indicated 
in column (14). The mean spectral index of each group is presented as well. * These stars are UX Orionis objects according to the definition of Dullemond et 
al. 2003. Note that the mean spectral index of these objects (3.13±0.03) is in good agreement with their believed group II nature (Dullemond et al. 2003). 
Measurement at 600 /um: 1.3±0.1 Jy for Elias 3-1 and 3.38±0.09 Jy for HD 163296. average of measurements by Ref. 15 and 17. ^ Based on the 1100 micron 
flux. " Measurement at 2 mm; ^ measurement at 2.6 mm: " measurement at 2.9 mm. 
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[350J 


[450J 


[750J 


[800J 


[850J 


[llOOJ 


[1300J 


12700J 




oJ 


n 




Ref. 


(1) 


K^) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




(11) 




(12) 


^ 1 Q\ 

(13) 


(14) 


V37d Cas 


I 


5.8±0.5 


1.3±0.2 


0.28±0.04 


0.23±0.02 




0.08±0.02 


0.038±0.006 


<0.0048 


8.8 


xlO" 


-'1 


4.59 


0.21 


3, 7, 13 


Elias 3-1^ 


I 


4.1±0.6 


2.6±0.4 


0.8±0.1 


0.78±0.05 


0.64±0.05 


0.35±0.05 


0.23±0.02 


0.05±0.01 


3.5 


xlO" 


-2 


3.22 


0.03 


1, 3, 7 


AB Aur 


I 


9±1 


3.8±0.6 


0.53±0.08 


0.50±0.05 


0.36±0.07 


0.15±0.02 


0.10±0.02 


0.011±0.005 


1.2 


xlO" 


-2 


4.32 


0.09 


3, 7, 8 


HD ,34282 


I 




1.6±0.2' 




0.41±0.03 


0.38±0.02 


0.18±0.02 


O.llitO.Ol 


0.024±0.003'' 


1.0 


xlO" 


-1 


3.36 


0.04 


9, 14, 15, 17 








0.22±0.04 






0.041±0.002 


0.04±0.01 


0.018±0.002 




1.2 


xlO" 


-2 


'X P-.A 

0.04: 




IK 17 1 Q 


















0.45±0.03 




9.2 


xlO'^ 


-2 


>2 94 






HD 100453 


I 














0.27±0.02 




2.1 


XlO" 


-2 


>2.82 




10 


HD 100546 


I 














0.47±0.02 


0.036±0.004'^ 


2.9 


xlO" 


-2 


4.20 


0.50 


5, 18 


HD 135344 


I 


6±2 


3.2±0.2 




0.57±0.02 


0.49±0.01 


0.21±0.01 


0.14±0.02 


<0.076'' 


1.6 


xlO" 


-2 


3.89 


0.05 


2, 15, 16 


HD 139614 


I 








0.61±0.03 




0.27±0.01 


0.24±0.01 


0.08±0.01" 


2.7 


xlO" 


-2 


3.17 


0.19 


15 


HD 142527 


I 








4.2±0.5 






1.19±0.03 


0.047±0.006'^ 


1.5 


xlO" 


1 


3.58 


0.41 


18, 19 


HD 169142 


I 




3.3±0.1 




0.55±0.03 


0.57±0.01 


0.29±0.01 


0.20±0.02 


0.07±0.02" 


2.4 


xlO" 


-2 


3.56 


0.11 


15, 16 


HD 179218 


I 














0.071±0.007 


0.008±0.001 


2.4 


xlO" 


-2 


4.06 


0.26 


9 


Mean 


I 






3.60 


0.09 




HD 31648 


II 




3.3±0.3 






0.78±0.02 




0.36±0.02 


0.04±0.01 


3.6 


xlO" 


-2 


3.44 


0.18 


8, 16 


UX Ori* 


II 














0.020±0.002 


0.0038±0.0004'' 


1.6 


xlO" 


-2 


3.20 


0.29 


12 


HD 35187 


II 




0.27±0.05 




0.12±0.02 


0.06±0.005 


0.08±0.01 


0.035±0.002 




5.0 


XlO" 


-3 


2.77 


0.38 


15, 17, 19 


BF Ori* 


II 














0.006±0.002 




1.0 


xlO" 


-2 


>2.82 




11 


HD 104237 


II 














0.09±0.02 




7.8 


XlO" 


-3 


>2.79 




5 


HD 141569 


II 




0.07±0.02 






0.014±0.002 


<0.036 


0.009±0.004 




5.6 


XlO" 


-4 


2.92 


0.38 


15, 17, 19 


HD 142666* 


II 




1.14±0.04 




0.35±0.02 


0.313±0.005 


0.18±0.01 


0.127±0.09 


<0.063'' 


1.6 


XlO" 


-2 


3.07 


0.01 


15, 16 


HD 144432 


II 








0.10±0.03 




0.07±0.01 


0.037±0.003 




4.9 


X 10" 


-3 


3.00 


0.62 


15, 19 


HD 150193 


II 




0.5±0.2 






0.12±0.02 




0.05±0.01 


0.014±0.003 


6.3 


XlO" 


3 


2.97 


0.11 


6, 8 


AK Sco* 


II 




0.24±0.03 






0.083±0.009 


0.036±0.004 






3.1 


XlO" 


3t 


3.05 


0.35 


6 


HD 163296'' 


II 


10.7±0.7 


5.6±0.4 


2.65±0.08 


2.32±0.04 


1.92±0.03 


1.02±0.05 


0.78±0.03 




6.8 


XlO" 


-2 


2.93 


0.09 


7 


WW Vul* 


II 














O.OlliO.OOl 




1.3 


XlO" 


-2 


>2.61 




11 


SV Cep* 


II 














0.007±0.002 




8.6 


XlO" 


-3 


>3.18 




11 


Mean 


II 






3.06 


0.06 
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the circumstellar disks of our sample sources based on 
the measured 1.3 mm flux, assuming the simple model of 
an optically thick disk with a black-body temperature of 
40K. The resulting radii are of the order of ~ 10-20 AU. 
Comparing this to the observed radii of Herbig Ae stars 
(e.g. Mannings & Sargent 1997, Natta et al. 2003, Testi 
et al. 2003, Pietu et al. 2003) or the typical disk sizes 
of T Tauri stars (~100 AU), the estimated radii of our 
sample-star disks seem too small. It is therefore unlikely 
that the disks of our sample stars are optically thick at 
(sub-)mm wavelengths. 

Supposing an optically thin disk at (sub-)mm wave- 
lengths and a typical grain temperature of 40K, one can 
estimate the disk mass from the 1.3 mm flux. Using for- 
mula (14) in Hillenbrand et al. (1992), we derived disk 
masses for the sample sources between 5.6 x 10~* and 

1.5 X 10~^ Mq. This compares well to the maximum disk 
mass of 0.3 M*, set by gravitational instabilities in the disk 
(e.g. Gammie 2001). When the disk is optically thick, its 
mass would be several factors higher, which in some cases 
(e.g. HD 34282) would be close to this theoretical limit. 

We have determined the value of the slope of the 
(sub-)mm SED by means of a power law fitted to the data 
between 350 and 2700 /im. It has not been proven that 
Herbig Ae/Be disks arc still optically thin at the shortest 
wavelengths of this interval. However, for the sources for 
which we have measurements at the shorter wavelengths, 
we can also compute the (sub-) mm slope based on the 
longer-wavelength data only. There is no significant dif- 
ference between the two. Hence, the slope of the energy 
distribution down to 350 and 450 micron does not deviate 
significantly from a power law. This shows that the optical 
depth at these wavelengths is not very different from that 
at 1.3 mm. Another way of saying this is that our results 
would not change when omitting the 350 and 450 micron 
data. 

Laboratory studies of the emission properties of dust 
species with astronomical relevance (e.g. Mennella et al. 
1998; Colangeli et al. 2001) show that variations in chem- 
ical composition only exert a very mild influence on 
the (sub-)mm emission characteristics of dust particles. 
Particle size, on the other hand, has quite a large influence 
on the (sub-)mm emission. Therefore, and in accordance 
with previous authors (e.g. Mannings & Emerson 1994; 
van den Ancker et al. 2000; Bouwman et al. 2000; Testi 
et al., 2003), we attribute the differences in (sub-)mm - 
spectral-energy-distribution slope to variations in dust 
particle size. 

It is clear from Fig. |21 and Table that the 
(sub-)mm spectral indices and the classification of the 
sample stars are correlated. Group I sources have an aver- 
age spectral index n equal to 3.60 ± 0.09. In group II, n is 

3.06 ± 0.06. Fig. 13 describes this correlation alternatively. 
The spectral index is plotted versus the flux excess at 60 
micron. The latter quantity is the difference between the 
measured flux and the flux of the reddened Kurucz model 
(see Fig. ^| at this wavelength. Group I objects, which 
are by definition brighter mid-IR sources, indeed display 
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Fig. 3. The (sub-)mm spectral index versus the flux ex- 
cess at 60 micron in magnitudes. Group I sources are rep- 
resented by filled, group II by open symbols. The upper 
dashed line represents the average value of n in group I, 
the lower line the average in group II. 



higher values for n. We interpret this correlation in terms 
of a correlation between the grain size distribution of the 
cold grains and the disk geometry. Group I sources (flared- 
disk systems) appear to have smaller grains than group II 
sources (self-shadowed-disk objects). 

Two explanations are given here. Due to coagulation, 
dust particles grow. It is therefore possible that different 
grain sizes indicate different evolutionary stages. Group 

I sources have dust that is only marginally evolved com- 
pared to ISM dust. Group II sources on the other hand 
have a significant population of mm-sized grains. We sug- 
gest that flared disks, in which coagulation takes place, 
evolve into geometrically somewhat flatter self-shadowed 
disks (by dust settling to the midplane) and become group 

II members. 

Alternatively, it is possible that, in contrast to coagula- 
tion, there are mechanisms that replenish the small-sized- 
grain population in circumstellar disks. The presence of 
the larger particles in self-shadowed-disk systems would 
then imply that this replenishment is not active/efficient 
(anymore) in these sources. In this hypothesis, flared 
disks are not necessarily evolutionally connected to self- 
shadowed disks, and the geometry of the disk is kept flared 
by the replenishment of the small grains. 

We note that the shape of the amorphous 10 micron 
silicate feature, which is linked to the particle sizes of 
the warm silicates (Bouwman et al. 2001, van Boekel et 
al. 2003), does not give any direct information about the 
grain size distribution of the silicates in the mid-plane of 
the circumstellar disk. This feature is believed to originate 
from the surface layers of the disk. Grain growth of these 
warm grains has been observed in this sample of HAEBEs, 
but is uncorrelated with the (sub-)mm spectral index (see 
AV04). Hence these small (~ 0.1-2 /im) warm grains seem 
to have an evolution, independent of the grain growth in 
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the interior of the disk. On the other hand, the 10 micron 
feature can be an indication for the sizes of the smallest 
grains in the disk, assuming that efficient vertical mixing 
in the disk occurs. 



lution grains will grow, this offers the first observational 
indications to the idea that flaring disks evolve into self- 
shadowed disks. 



4. Conclusions 
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We showed in this paper that there is a significant dif- 
ference between the (sub-)mm spectral indices of HAEBE 
stars with different IR SEDs. This could very well mean 
that the grain size distribution of the cold grains in a 
HAEBE star depends on the disk geometry. 

The HAEBE systems with large IR excesses (group I, 
i.e. flaring-disk objects) have smaller cold grains than sys- 
tems with more moderate IR excesses (group II, i.e. self- 
shadowed-disk sources). In other words, it appears that 
the flaring geometry of a circumstellar disk can only per- 
sist when a sufficiently large population of small grains 
exists in the disk. 

The circumstellar-disk models of DD04 agree qualita- 
tively with the results presented here. Disks in which the 
vast majority of the grains have grown to larger sizes, lead 
automatically to self-shadowed geometries. Models with 
small grains result in a mid-infrared SED akin to group I 
sources. Models with grain sizes larger than ~1 mm not 
only display mid-infrared SEDs more resembling those of 
group II sources, but also change their (sub-)mm slope 
from larger values of n to a black-body character (n = 3). 
DD04 note that the (sub-)mm slope changes with some- 
what smaller grain sizes than are required to alter the 
mid-IR SED from group I to group II. Therefore, accord- 
ing to the models, the spectral index can have all values 
larger than 3 in group I sources, and all group II sources 
should have a black-body (sub-)mm slope. This is con- 
sistent with the observed trend in which group I sources 
are observed to have 3.17 < n < 4.59 and group II sources 
have n very close to 3. The DD04 models predict that grain 
growth resulting in a change of n from values similar to 
those found in the ISM (> 3) to the optically thin black- 
body value (n = 3) occurs before the disk becomes self- 
shadowed. However the DD04 models do not model the 
dynamical processes leading to grain growth. At present 
the time scale of this transition from larger n to n = 3 
is unknown, so we are unable to estimate the fraction of 
group I sources that should have n = 3. Therefore we 
have no reason to assume a priori that our distribution of 
(sub-)mm-spectral-index values is incompatible with the 
interpretation outlined in Section 13.21 

In the present paper, we report a clear correla- 
tion between the shape of the mid-IR SED and the 
(sub-)mm slope in our sample. Previous authors have 
shown that variations in (sub-)mm slope within systems 
with proto-planetary disks are caused by variations in 
dust particle size. D02 have shown that the shape of the 
mid-IR SED is a proxy for disk geometry (flared versus 
self-shadowed). In accordance with these results, we inter- 
pret our find as a correlation between disk geometry and 
grain size. If one assumes that in the course of disk evo- 
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